Absence of decoherence in the complex potential approach to nuclear scattering 
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Time-dependent density-matrix propagation is used to demonstrate, in a schematic model of an 
open quantum system, that the complex potential approach and the Lindblad dissipative dynamics 
are not equivalent. While the former preserves coherence, it is destroyed in the Lindblad dissipative 
dynamics. Quantum decoherence is the key aspect that makes the difference between the two 
approaches, indicating that the complex potential model is inadequate for a consistent description 
of open quantum-system dynamics. It is suggested that quantum decoherence should always be 
explicitly included when modelling low-energy nuclear collision dynamics within a truncated model 
space of reaction channels. 
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Quantum coherence manifests itself through quantum 
interference effects in quantum dynamical systems, which 
is clearly demonstrated in double-slit experiments with 
single electrons. This property of matter is diminished 
(decoherence) when an external environment interacts 
with the system [i|-y]- Quantifying the role and impor- 
tance of decoherence in quantum systems is now perva- 
sive in various branches of physics and chemistry, includ- 
ing studies of quantum measurement and quantum in- 
formation [3] ■ The concept of a reduced (but not closed) 
quantum system evolving in the presence of couplings to 
an environment of complex states is common throughout 
disciplines. 



The schematic model refers to the scattering of a wave- 
packet Tpkoi^) (with the mean wave- number fco) off a po- 
tential barrier ^(a;). This is considered a reduced quan- 
tum system, which irreversibly interacts with a complex 
environment of excluded degrees of freedom. Of inter- 
est is the effect of such an environment on the reduced 
quantum-system dynamics. It is here investigated us- 
ingtwo descriptions: (i) the complex potential approach 
[3-[3] , and (ii) the Lindblad dissipative dynamics [^, Q . 
The key question addressed in this paper is: Are these 
approaches equivalent? A common view in the nuclear 
physics community is that they are. This general belief 
inhibits research into quantum decoherence which is at 
the same time highly topical in other areas of physics 
and chemistry. My model calculations clearly demon- 
strate that quantum decoherence makes the difference 
between the two descriptions, and is not accounted for 
in the complex potential model. The methodology will 
be highlighted first. Afterwards, the calculations are pre- 
sented and discussed, and finally a summary is given. 

The Liouville-von Neumann master equation dp/ dt = 
C p dictates the time evolution of the reduced-system 
density-matrix operator p{t). Its initial value describes a 
pure state, and is determined by the initial wave-packet 
as yo(0) = |'0)('!/'|. The Liouvillian C p \s different in the 
two descriptions studied: 



(i) In the complex potential approach j^], this is 



Cp 



.HeffP - pHLf), 



(1) 



where the effective non-Hermitian Hamiltonian 
Heff = Hs-iW{x), being H^ = f + V{x) the Her- 
mitian Hamiltonian of the reduced system (T is the 
kinetic energy operator) and W{x) > describes 
the irreversible environmental interaction(s). 

(ii) Irreversibility in dynamics of an open quantum sys- 
tem can be consistently described by the Lindblad 
master equation [8] . Here the Liouvillian reads as 



Cp = - 



H,, 



E(C'"'5Cl-J[ClCa,p]J 



(2) 



where [. . .] and [...]+ denote the commutator and 
anti-commutator, respectively. Each Ca is a Lind- 
blad operator for a dissipative coupling, physically 
motivated according to the specific problem. We 
use a spontaneous emission Lindblad operator jlO|, 
^21 — •\/7li|2)(l| that describes a decay from 



I.e. 

the reduced-system state |1) to the environmental 
state 1 2). (These states are assumed to be orthonor- 
mal.) State |2) is not a reaction channel, but an 
auxiliary state [10| supplying a probability drain 
only. This state mocks up a high density of com- 
plex states (environment) and irreversibly absorbs 
probability from the reduced system. The absorp- 
tion rate to state |2) is given by 7^;|. = W{x)/h 
where W{x) is taken as in approach (i). Thus, ap- 
proaches (i) and (ii) refer to a single reaction chan- 
nel problem. 

The Liouville-von Neumann master equation with ei- 
ther III) or © is then represented in a grid-basis [ll| . The 
density-matrix elements obey the following equations: 

(a) In the approach (i), 
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where 



{Cdp)xx' = -j^{W{x)+W{x'))p^ 



(b) Whereas in the approach (h), 
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(5) 
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= (Cdp) 
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and the elements 



and 



(6) 
obey equations of 



motion hke ([6]). The dissipative terms are 
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where (kj = 1,2) and 7L' = E.^klix M- K 
makes the second term of ([7]) equal to (|3]) for the 
elements p]}^, that describe the reduced system. 
Thus, the first term of (O makes the difference 
between the two approaches, which is essential as 
shown below. It is worth mentioning that the el- 
ements (0^^, absorb probabilities only, whilst pl^., 



and 



rxx' 

\ikepll,. 



are always zero. These are initially zero 



Having obtained the dynamical evolution of the density 
matrix, expectation values of an observable O result from 
the trace relation (0'(i)) = Tr((!)/5(i)). 

Results and discussion. In the model calculations the 
grid (-50 to 50 fm) was evenly spaced with 256 grid- 
points. The potentials V{x) and W{x) are the Eckart 
potential jl3| and a Gaussian function, both centered 
at zero, i.e., V{x) = Vq cosh"'^ {x / ao) and W{x) — 
Wo exp[—x-^ /2aQ^]. The parameters are (Vo,ao) = (61.1 
MeV, 1.2 fm) and {Wo,aow) = (5 MeV, 1.5 fm). A 
minimal-uncertainty Gaussian wave-packet describes the 
relative motion of two objects with mass numbers Ai = 
16 and A2 = 144. The wave-packet is initially centered at 
xq = 25 fm, with width ao = 5 fm, and was boosted to- 
wards the potentials with the appropriate average kinetic 
energy for the initial total energy Eq = 55 MeV. 

The density-matrix was propagated in time using the 
Faber polynomial expansion of the time-evolution oper- 
ator |14|, and the Fburier method 11] for the commu- 
tator between the kinetic energy and density operator. 
The time step for the density-matrix propagation was 
Ai = 10^^^ s. The numerical accuracy of the propa- 
gation was checked using a time-dependent calculation 
without the dissipative term Cd p in ^ 01 ^. It was 
confirmed that the trace and purity of the density-matrix, 
Tr(p) = Tr(/5^) = 1, and the expectation value of the 
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FIG. 1: (a) The energy mean value of the reduced system as a 
function of time in the complex potential approach (solid line) 
and in the Lindblad dissipative dynamics (dashed line), (b) 
The same but for the trace of the density matrix. Lindblad's 
dynamics preserves the total probability (dotted line), i.e., 
the sum of the probability in the reduced system and in the 
environment. 



system energy Tr (Hp) were maintained with high accu- 
racy over the required number of time steps, typically 130 
when the centroid of the recoiled body of the wave-packet 
reaches 25 fm. 

Figure [T] shows the time evolution of the reduced- 
system energy (Panel a) and the trace of the density- 
matrix (Panel b) using both the complex potential ap- 
proach (solid line) and the Lindblad dynamics (dashed 
line). These quantities show similar features, their val- 
ues declining in time as expected. The total probabil- 
ity is mantained in Lindblad's dynamics (dotted line), 
whilst the complex potential approach cannot guaran- 
tee its preservation. The disagreement of the solid and 
dashed lines is due to the first term of ([7]). Removing it, 
the two approaches provide the same results. However, 
this term is crucial, making the two approaches physically 
different due to the quantum decoherence, as shown in 
Fig.H 

The measure of coherence in the reduced system is the 
ratio Tr(/92)/[Tr(/5)]2 [Hill]. Its time evolution is pre- 
sented in Fig. [2]for the complex potential approach (solid 
line) and the Lindblad dynamics (dashed line). The for- 
mer preserves coherence, whilst the Lindblad dynamics 



results in decoherence. The complex potential leads to 
a coherent damping of all the elements of the density 
matrix, affecting only the amplitude of the off-diagonal 
elements. However, the first term of ([7]) results in dephas- 
ing, destroying the initial phase relationship between the 
off-diagonal elements. This is quantum decoherence. 

An interesting question arises: Does decoherence affect 
observables such as the quantum tunnelling probability? 
Yes, it does, as shown in Fig. [31 Here, the tunnelling 
probability through the Eckart potential barrier is pre- 
sented as a function of the initial total energy Eq. This 
probability is numerically calculated (after a long propa- 
gation time) as the trace of the density-matrix over neg- 
ative values of the x coordinate, i.e., Tr(P/5) where the 
projector P = 1 — ©(a^), being Q{x) the Heaviside step 
function. The dotted line represents the analytic trans- 
mission coefhcient [13], whilst the thin solid line shows 
the tunnelling probability for the wave-packet with a def- 
inite energy spread. This results in the small difference 
between the two lines. With respect to these, the tunnel- 
ing probability is reduced in both the complex potential 
approach (thick solid line) and the Lindblad dynamics 
(dashed line). The difference between these two lines 
demonstrates that decoherence substantially changes the 
quantum tunnelling probability, i.e., up to a factor ten at 
the lowest incident energies. The energy dependence of 
this probability is clearly affected by decoherence (com- 
paring the thick solid and dashed lines with the thin 
solid line). These results provide quantitative support 
to the idea [l7Hl9| that explicit inclusion of quantum de- 
coherence may well be necessary to consistently describe 
(within a truncated model space) low-energy nuclear col- 
lision dynamics. 

Quantum decoherence is a key aspect of irreversibil- 
ity, which is mostly overlooked in quantal models of low- 
energy nuclear phenomena. It should be treated simul- 
taneously with energy dissipation and coherent quantum 
superpositions, the latter being the basis of the coherent 
coupled-channels approach to nuclear reaction dynamics 
around the Coulomb barrier. Coherent quantum super- 
positions manifest themselves through the measurement 
of the experimental fusion barrier distributions, while 
energy dissipation is revealed in deep-inelastic scatter- 
ing that also occurs at near-barrier energies fl9l . While 
the coherent coupled channels calculations [20| are able 
to explain several collision observables, major problems 
are unresolved. Foremost is the inability to describe the 
elastic and quasi-elastic scattering and fusion processes 
simultaneously [2l| and the related, more recent fail- 
ure to describe consistently below-barrier quantum tun- 
nelling and above-barrier fusion yields [i7|. New, pre- 
cise measurements have inevitably led to phenomenolog- 
ical adjustments [23, |23| (sometimes contradictory) to 
stationary-state coupled channels models to fit the exper- 
imental data, but without a physically consistent founda- 
tion. The present work suggests that quantum decoher- 
ence should be taken into account seriously in modelling 
fusion and nuclear scattering. This may also have signifi- 
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FIG. 2: The same as in Fig. \T\ but for the measure of coher- 
ence [I^l in the reduced system. While the complex potential 
approach preserves coherence (solid line), it is destroyed in 
the Lindblad dissipative dynamics (dashed line) . 
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FIG. 3: Tunnelling probability through the Eckart potential 
barrier as a function of the initial total energy: analytic prob- 
abilities (dotted line), probabilities of the wave-packet with a 
definite energy spread (thin solid line), and the latter within 
the complex potential approach (thick solid line) and the 
Lindblad dissipative dynamics (dashed line). The environ- 
ment reduces the tunnelling probability, playing decoherence 
a crucial role (comparing the thick solid line with the dashed 
line, and these with the thin solid line). 



cant implications for open systems-related approaches to 
phenomena in near-threshold exotic nuclei, such as con- 
tinuum shell models [24] and approaches to understand- 
ing the quenching of spectroscopic factors [25l . |26| . 

In summary, the non-equivalence between the com- 
plex potential approach and the Lindblad dissipative dy- 
namics for describing the time evolution of a schematic, 
open quantum system has been demonstrated. Time- 
dependent density-matrix propagation shows that quan- 
tum decoherence makes the difference, which is not 
treated in the complex potential model. It is inadequate 
for a consistent description of open quantum system dy- 
namics. Decoherence substantially affects the quantum 
tunnelling probability, and should be included in a com- 



plete description of low-energy nuclear scattering. 
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